Extended Debye Model for Molecular Magnets 
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Heat capacity data on Mni2 are fitted within the extended Debye model that takes into account 
a continuum of optical modes as well as three different speeds of sound. 
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Molecular magnets (MM) such as Mn.12 (Ref. 0) are 
relatively new materials that have a giant effective molec- 
ular spin S (such as S = 10 for Mni^) built from several 
atomic spins by a strong intramolecular exchange interac- 
tion. The magnetic molecules have a uniaxial anisotropy 
that is responsible for magnetic bistability and long re- 
laxation over the barrier at low temperatures^ As the 
magnetic core of these molecules is surrounded by or- 
ganic ligands, the exchange interaction between different 
molecules building a crystal lattice is very small. This 
allows them to relax independently from each other, in 
contrast to ferromagnets. A fascinating phenomenon dis- 
covered in molecular magnets is resonance spin tunneling 
under the barrier that happens if the energy levels of the 
spin S in both wells matchi 4 ^^ 

Molecular magnets is a new type of condensed mag- 
netic systems whose properties differ from those of ferro- 
magnets and dilute paramagnets. Although in the most 
temperature range MM are paramagnetic, their relax- 
ation can differ from that of a single spin embedded in 
an elastic matrix. Since the wave length of emitted and 
absorbed phonons or photons exceeds the lattice spacing, 
there can be pronounced coherence effects in relaxation 
such as superradiance.— Photon^ and phonon- superra- 
diance in MM can increase relaxation rates by a huge 
factor. On the other hand, the opposite effect for ini- 
tial states of spins with random phases should lead to 
suppression of the rates by a huge factor. Strong inho- 
mogeneous broadening in MM tends to destroy coher- 
ence effects, however, so that efforts should be done to 
understand the relaxation data. Another collective phe- 
nomenon in relaxation that is not yet fully understood 
theoretically is the phonon bottleneck (see Refs. [TolflTI for 
older reference and Refs. [l2lfl3l for recent work). 

To be able to test more sophisticated collective models 
of relaxation in MM, one should have reliable theoretical 
estimations of the single-spin relaxation rates, most no- 
tably the one-phonon or direct relaxation rate. The latter 
can depend, in general, on spin-phonon couplings that are 
difficult to measure. On the other hand, there is a simple 
mechanism of spin-lattice coupling through rotations of 
the magnetic molecules by transverse phonon o 14 : 15 that 
can serve at least as the low bound on spin-lattice relax- 
ation. As in this mechanism the crystal field acting on 
the spin is not distorted but only rotated, no unknown 
coupling constants enter the theory. Also this mechanism 
is likely to be the dominating relaxation channel since 
the cores of magnetic molecules should be much less de- 



formable than the ligands. The corresponding results for 
the relaxation rates T due to direct processes, as well as 
the Raman processes, depend on only one parameter that 
is currently not precisely known, the speed of transverse 
sound Vt- For direct processes one has T oc 1/uf , whereas 
for Raman processes T oc l/v\° . Thus the uncertainties 
of Vt dramatically amplify in the relaxation rates. 

In the absence of direct measurements of the speed of 
sound, the latter can be extracted from the heat capacity 
of the lattice by fitting the measured C(T) to the Debye 
theory and extracting the Debye temperature Od that is 
proportional to the speed of sound. So Od = 38 K of Ref. 



16 results in v = 1600 m/s, whereas Od = 41 K of Ref. 
171 results in v = 1727 m/s. In fact, determination of Od 
in Refs. [I6IIT7I relies on the low-temperature data, where 
the phonon contribution to the heat capacity is C p h oc T 3 
and the Debye model with the rigid the cut-off at the 
Brillouin-zone boundary, that is a crude approximation, 
is not actually used. 

A problem with the extracting the Debye temperature 
and the speed of sound above is the assumption that the 
three branches of acoustic phonons in the crystal have 
the same speed v. Indeed, at low temperatures one has 
Cph oc (vi 3 + i'2 3 + w^ 3 ) T 3 , and there is no way to find 
Vi separately from this formula. On the other hand, 
acoustic phonon modes in such a complicated crystal 
as Mni2 should not be strictly longitudinal and strictly 
transverse and all three speeds of sound, as well as all 
three Debye temperatures, should be different. If they 
all differ much, the smallest of them dominates the low- 
temperature heat capacity and, to a much greater extent, 
the spin-phonon relaxation rates. If one assumes that 
there are two degenerate transverse phonon modes with 
v% = v% = Vt <C U3 = Vi, then instead of Od = 38 K and 
v = 1977 m/s one obtains ©o,t = 33 K and v t — 1717 
m/s. If there is one phonon mode with v± <C «2,3> then 
one obtains ©d,i = 26.3 K and v\ = 1368 m/s. 

To improve the description and extract more infor- 
mation, one can use the heat-capacity data at higher 
temperatures where the contributions from the phonon 
modes add up in a different way. This is where the Debye 
model begins to really work and crudeness is introduced 
into the theory. In addition, optical modes become very 
important and should be taken into account. In recent 
Ref. [H on another MM Fes the Debye model has been ex- 
tended by adding an Einstein oscillator that accounts for 
all optical modes, and the analysis in a broader temper- 
ature range rendered Od = 19 K (as well as the Einstein 
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FIG. 1: Heat capacity of Mn.12. The standard Debye model 
fails everywhere except for very low temperatures, while the 
extended Debye model perfectly fits the data. 



temperature B = 38 K) that differs much from the pre- 
viousely extracted value ©d = 34 K for this MM. On the 
other hand, all three speeds of sound were considered as 
the same in this a naly sis. 

The data of Ref. Ha shown in Fig. Q] go as C oc T at el- 
evated temperatures, and the heat capacity per molecule 
by far exceeds the value 3&b of a crystal with one atom 
per unit cell at T > ©d- This means that there are a 
lot of optical modes forming a continuum with a nearly 
constant density of states. This is expected for molec- 
ular magnets having hundreds of atoms within the unit 
cell. Description based on a continuum of optical modes 
is much more reasonable than the Einstein theory with 
all due respect. The aim of the present paper is thus to 
formulate the extended Debye model (EDM) including a 
continuum of optical modes. This will be used to extract 
the speeds of sound in Mni2 from the experimental data 
without making the assumption v% — vi = V3. 

The thermal energy of acoustic phonons per unit cell 
of a crystal is given by 
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where N is the number of unit cells in the crystal, WkA are 
phonon frequencies and A = 1,2,3 are phonon polariza- 
tions. One can replace summation over k by integration. 
Within the Debye model one assumes that the relation 
k>kA = v\k holds everywhere in the Brillouin zone that 
is approximated by a sphere bound by the Debye wave 
vector fcrj ■ The latter is defined by the requirement that 
that total number of phonon modes is N : 
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where vq is the unit-cell volume. This yields 
kn = (67r 2 /w ) 1/3 . 



(3) 



One can introduce Debye frequencies and Debye 

temperatures Od.a for different acoustic phonon branches 



A as 

^D,A = V\kr>, fc B 6D,A = MIb,X- (4) 

Now Eq. ([1]) can be written as 
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where the densities of states are given by 

p x (w) = 3w 2 /Q^ x (6) 

for to < f2rj,A and zero otherwise. At low temperatures, 
T<0 D , Eq. © yields 
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At high temperatures, T ^> ©d, Eq. §5§ yields 

U ph = 3k B T, C ph = 3fc B - (9) 

For all 0d,a being the same, the coefficient in Eq. ([8|) is 
a huge number 127r 4 /5 ~ 234. Because of this, Eq. ([8|) 
does not smoothly join with Eq. © at T ~ 0d, and 
the applicability of Eqs. ([7]) and (SJ requires in fact very 
low temperatures, not just T <C ©i> On the other hand, 
Eq. (J9j> is at striking contradiction with the experiments 
shown in Fig. [1] because of the huge unaccounted con- 
tribution of optical modes. Thus the usefulness of the 
Debye model in its standard form is limited, at least for 
molecular magnets. 

To improve the Debye model in a minimal way, one 
can add optical modes with a constant density of states 
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where J7g is another characteristic frequency that should 
be considered as a fitting parameter. With Eqs. ([6]) and 
(fT0|) inserted, Eq. (© becomes 
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where /cb@g = fi^G and 
F n (y) = 
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For the heat capacity one obtains 
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FIG. 2: Heat capacity of Mni 2 in C(T) /T vs T 2 representa- 
tion used to eliminate the small parasite term C(T) cc T. 
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FIG. 3: Heat capacity of Mni2 in the presentation used as 
a crucial test for different models. One can see that only 
the model with three different acoustic phonon modes fits the 
data. 



In the high-temperature limit T 3> 8d these equations 
yield 

C/ph = C ph - 37r 2 fc B? ^, (14) 

instead of Eq. @ and in accord with the experimental 
data of Ref. [l6| shown in Fig. [TJ 

To fit the experimental heat capacity with Eq. (|13l) . 
one has first to subtract the spin (Schottky) contribution 
Cs{T) from C(T). Using the spin Hamiltonian in zero 
field 

H = -DS 2 z - ASt + H', (15) 

where D/k B = 0.548 K, A/k B = 1.1 x 10~ 3 Kj±&22 and 
H' is the part of the Hamiltonian that does not com- 
mute with S z , nonessential for Cg(T) in Mni2. With the 
spin levels e m — —Dm 2 — Am the energy is given by 
U s = (l/Z)J2 S m= „ s e m exp[-e m /(k B T)}, where Z = 

E™ = -s ex P h«W (k B T)] . Then C s = dU s /dT, and fit- 
ting the lattice heat capacity C p h = C — Cs yields the 



EDM curve shown in Fig. [T] that is in an excellent ac- 
cord with the experimental data in the entire tempera- 
ture range. 

In fact, the lattice heat capacity also contains the con- 
tribution of nuclear spins ~ 1/T that is small in the 
Kelvin range, as well as the T-linear term that must be 
an artefact of the experimental procedure. To eliminate 
the parasite T-term and visualize the C p h oc T 3 depen- 
dence at low temperatures, it is convenient to plot C p h/T 
vs T 2 to get a straight line, as shown in Fig. [5] The slope 
of the straight line yields the average value of the Debye 

temperature 6d = ©d^/S^ = 38 K, in accord 
with Ref. UH, whereas the coefficient in the T-linear term 
is 0.022. Subtracting this small constant and fitting the 
rest, one obtains 

6 D ,i = 29.6 K, e D:2 = 47.8 K, 9 D:3 = 61.0 K, (16) 

as well as 

9 G = 15.6K. (17) 

From Eqs. ((4]) and (fTB]) follows 

Vi = 1541 m/s, v 2 = 2488 m/s, v 3 = 3176 m/s. (18) 

The first of these speeds of sound should correspond to a 
nearly transverse mode, while the last one should corre- 
spond to a nearly longitudinal mode. The former is the 
most important in relaxation. 

Of course, one can say that with enough fitting param- 
eters one can fit any function. While in general it is true, 
the scheme used here is a minimal model with no exces- 
sive fitting parameters. Accounting for optical phonons 
with a single parameter Qq is a must, and there are no 
physical reasons to set speeds of sound the same. The ex- 
perimental results in the natural representation in Fig. [1] 
can also be fitted by the EDM with one acoustic phonon 
mode (vi — v% — v 3 = v) and two acoustic phonon modes 
(vi = v% = Vt, U3 = vi), and the results are visually not 
dramatically worse. However, this fitting method is in- 
ferior since it puts more weight on the high-temperature 
range and tends to neglect the low-temperature range. It 
is better to fit the C/T data with the subtraction of the 
parasite T-term, as was done above (see Fig. [5]). Then 
one can see more difference between different fits. The 
most stringent check of different fits can be achieved in 
the most balanced representation of C(T) over the whole 
temperature range using C/T 2 (with subtracted T-term) 
as the fitting target. The results shown in Fig. [3] demon- 
strate that only the model with three different speeds of 
sound really fits the data. In addition, for this model re- 
sults obtained with different fitting methods do not differ 
much. Fitting of C/T 2 with three phonon modes yields 
6 D ,A = 29.7, 47.1, 61.4 K and 6 G = 15.7 K that is very 
close to the results of the C/T fitting, Eqs. $TB]) and fI7). 
and even to the results of the C fitting (0d,a = 29.5, 50.0, 
57.0 K and Qq = 15.8 K). To the contrary, models with 
one or two different speeds of sound yield very differ- 
ent results with the three different fitting schemes. Thus 
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FIG. 4: Phonon density of states for M1112 within the ex- 
tended Debye model. 



one concludes that these models do not work. Still we 
quote the results from fitting C/T 2 in Fig.[3]for the refer- 
ence: One phonon mode 9d = 39.4 K and 9q = 18.4 K, 
two phonon modes 8d,* = 38.1 K, 9d,! = 41.3 K, and 
6 G = 17.5 K. 

Fig. |4] shows the density of states in Mni2 within the 
extended Debye model with parameters given by Eqs. 
dHJ) and (17). Although there are steps in the DOS 
that reflect the crudeness of the underlying Debye model, 
there are three smaller steps instead of a single large step 
in the original Debye model. Thus the EDM is much 
more realistic than the DM. The accuracy of the ex- 
tracted data on speeds of sound is difficult to estimate 
because of the assumption of the rigid cut-off at &d- One 
rather should consider Eqs. (fTrJ|) - (fT5)) as qualitative re- 
sults that capture some physics of phonons in MM. 

A practical question is how to apply the results ob- 



tained above to the relaxation in MM. All existing formu- 
las are based on the model with one longitudinal phonon 
mode and two transverse phonon modes. Within the 
mechanism of the molecule rotation without distortion, 
the contributions of the processes contain the factors 
sin 2 6 = [k x ekA] 2 /k 2 , where ekA are phonon polariza- 
tion vectors. Thus longitudinal phonons do not make a 
contribution while transverse phonons make the maximal 
contribution. In reality phonons are not purely longitu- 
dinal and purely transverse and it is difficult to extract 
the angle between k and ekA- For an estimation one can 
propose a rule of thumb: 8 = 7r/2, ir/4,0 for the phonon 
modes with vi, v%, and V3, respectively. With this con- 
jecture, one can make the following replacement in the 
formulas for the rates of direct processes: 
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For the values listed in Eq. (jl8|) . the correction term in 
square brackets is only 0.046 and can be neglected. Thus 
the rule of thumb simplifies to keeping only the softest 
mode with the coefficient 1/2. Then the increase of the 
rate due to using the model with three different phonon 
modes instead of the traditional model with one longi- 
tudinal and two degenerate transverse phonon modes is 
given by wf/(2u?) = 6^/(20^). With 6 D ,i = 29.7 
and Oo,t = 38.1 K obtained above, the rate increase 
makes up 6^/(29^) = 2.4. 
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